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Abstract - Motivated by the spin self-rephasing recently observed in an atomic clock, we intro- 
duce a simple dynamical model to study the competition between dephasing and synchronization. 
Two spins S are taken to be initially parallel and in the plane perpendicular to an inhomogeneous 
magnetic field A that tends to dephase them. In addition, the spins are coupled by exchange 
interaction J that tries to keep them locked. The analytical solution of the classical dynamics 
shows that, there is a phase transition to a synchronized regime for sufficiently large exchange in- 
teraction J > A compared to the inhomogeneity. The quantum dynamics is solved analytically in 
four limits - large/small J/ A and large/small S - and numerically in between. In sharp contrast 
to the classical case, the quantum solution features very rich ^-dependent multiscale dynamics. 
For any finite 6, there is no synchronization but a crossover around J = A between two regimes. 
The synchronization transition is only recovered when S oo, approaching the classical solution 
in a non-trivial way. 
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^ Introduction. — Synchronization is a collective phe- 
^nomenum that occurs in the dynamics of many different 
systems, see e.g. |1 and references therein. A famous 
if) classical model of synchronization is that introduced by 
^-^Kuramoto [2j|3]. In this model, two or more oscillators 
^^with distinct frequencies can synchronize when they are 
CO coupled by a sufficiently strong non-linearity. Recently, 
H a synchronization transition was observed in a quantum 
systems made of a large ensemble of spins 1/2 [4j|5] . In an 
V"! atomic clock, N ~ 4.10 4 trapped two-level atoms behav- 
rNing as pseudo-spins were found to synchronize beyond a 
critical density or interaction strength. In the experiment, 
the contrast of the Ramsey fringes measures the coherence 
of the atomic clock and typically decays in time. It was 
found that this decay almost stops at sufficiently large den- 
sity, substantially increasing the coherence of the clock. A 
simple picture explaining this self-rephasing was proposed 
in El. It involves two equal populations of atoms - cor- 
responding to hot and cold atoms - each represented by 
a macrospin S = N/4. The two macrospins feel different 
longitudinal magnetic fields, because of the spatial inho- 
mogeneity of the atomic cloud. In addition, atom- atom 
collisions generate an effective exchange coupling for the 
two macrospins. When the latter is strong enough, it im- 



pedes dephasing of the macrospins. 

Building further on this picture, we introduce a quan- 
tum model of two spins 6 to describe the dynamical com- 
petition between dephasing by the inhomogeneous mag- 
netic field and synchronization by exchange interaction. A 
legitimate question to ask is whether a genuine synchro- 
nization transition exists in a quantum model 6 . Here, 
we answer that question by presenting our results, leaving 
details of derivations to a longer companion paper [7l. We 
start by introducing the model, before giving its classical 
solution and several approximate quantum solutions. Our 
aim is to provide a qualitative picture of the quantum dy- 
namics as a function of the spin size S and of the ratio 
between exchange and inhomogeneity. 

Two spins model. — We consider the dynamics of 
two spins 61,62, of size 6, coupled by exchange interac- 
tion and subjected to an inhomogeneous magnetic field in 
the z direction. The corresponding Hamiltonian ressem- 
bles that of the two- level BCS model [8] and reads H = 
J s 6i-62+A s (6f — 6|) with J s (resp. A s ) the characteristic 
exchange (resp. inhomogeneity) energy (hereafter h = 1). 
The initial state is taken to be a coherent state in the direc- 
tion x perpendicular to the magnetic field: 1^(0)) = \Sf = 
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5, 5f = S). We characterize the state \i/j(t)) at time t by 
three quantities: the single spin contrast C\(t) = ^ Sl ^ 2 ^ ; 
the single spin unit vector direction H\ 2(f) = ,)J" 2 {, 

and the total spin contrast C(t) = ^ Sl ^ S ' 2 ^ , the lat- 
ter corresponding to the quantity that is measured in 
Ramsey fringes exper iments [ 4|[5] . At any time we can 
write C(t) = Ci(t)y/ (1 + Hi ^2)/ 2. For classical spins 
C\{t) = 1 since Si^(t) = SHi^(t). For quantum spins 
C\ (t) is in fact a direct measure of the effective spreading 

width Di(t) = ^ Sl ~^2 1 ^ - of the single spin wavepack- 
ets; indeed, quite generally C\(t) = 1 — (D\(t) — §)• For 
our specific choice of initial state, the initial values are 
C(0) = Ci(0) = 1, Di(0) = § and H ly2 (0) = e x ; further- 
more at any time we also have H 2 (t) = [nf , — nf, — nf] such 
that we can rewrite C(t) = Ci(£)|nf 
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Fig. 1: Contrast C as a function of rescaled time t [units of ^] 
for S = 10. Comparison between classical dynamics (dashed 
green line) , quantum dynamics (dotted blue line) and cumulant 
approach (full red line) in 3 regimes: (a) J/A = 0.2 (dephased); 
(b) J/A = 1 (critical); (c) J/A = 10 (synchronized). Initial 
decay (U), fast oscillations (t ) and envelope (t e ) times are 
indicated. 

Classical synchronization transition. — We first 
consider the classical version of the two spins model. We 
get rid of the spin size S by defining rescaled quanti- 



ties: time t = A s 1 t, exchange J = A 2 J s , inhomogeneity 
A = A S A S and spins Hi^ 2 = A^Si^ with A s = S. Writ- 
ing Hi = [cos#, sin#^l — (^) 2 sin 2 ^ sin 2 0], the spin 
equations of motion ni = [Ae 2 + Jn2(t)] x ni(f) map to a 
nonlinear pendulum equation 0(t) = A^/l — ^| sin 2 

[7]. The solution #(£) is obtained as the incomplete ellip- 
tic integral of the first kind F(sm6; £) = At such that 
the contrast is C(t) = |cn[At; ^]| with cn the Jacobi el- 
liptic function. At short time, 6(1) ~ At is linear in time 
such that C(i) ~ 1 - ^ w 1 - P/ff with k = For 
larger time the dynamics is nonlinear and, depending on 
the ratio ^, we distinguish three regimes. The case J < A 
corresponds to the dephased regime; the contrast vanishes 
periodically with a period t p = ^K(^) « J with if(fc) 
the complete elliptic integral of the first kind (Fig. [T^l, 
green dashed). The case J > A corresponds to the syn- 
chronized regime, which, by definition, means that the 
contrast remains finit e at all times. It reaches its min- 
imal value C m i n = ^Jl — ^ periodically with a period 
i P = 7#(y) ~ 7 (Fig. [j]:, green dashed). For small 
| J — A | the period t p diverges as t p ~ ^ log rj~A| ~~ ^ 00 
signaling a critical regime at J = A. At this point the 
contrast C(t) = sech(At) monotonically decreases on the 
short timescale U = ^ (Fig. [IJd, green dashed). 

To complete the characterization of the classical spin 
dynamics, the trajectory of the single spin direction 
H\(t) is plotted in Fig. [3] (full red line) for each 
regime. For later comparison with quantum results, 
we quote the analytical asymptotic results: H\(t) « 
[cos At , sin At , 0] when ^ <C 1, and H\(t) « [1 — 
$s sin 2 (Jt), f sin(Jt) cos(Jt), f sin 2 (Jt)] when ^ > 1. 

Below, when discussing quantum dynamics, we shall 
still refer to the three regimes as dephased (J « A), crit- 
ical ( J ~ A) and synchronized (J ^> A), even if synchro- 
nization stricto sensu does not occur. 

Quantum dynamics. — We now consider the quan- 
tum dynamics. In order to conveniently compare with the 
classical results, we also define the rescaled quantities t, J 
and A, however now A s = \/S(S + 1). [^] We emphasize 
that after rescaling the quantum model still explicitly de- 
pends on the spin size S since the latter fixes the size of the 
Hilbert space as (25 + l) 2 . In other words, the quantum 
dynamics is controled by two dimensionless parameters: 
J/A and S. In the following, we study the quantum dy- 
namics of the three quantities of interest C(t), Ci(f), H\(t) 
using five different and complementary approaches. 

Before going to the quantitative results, we first sketch 
the typical behavior of the quantum contrast C(t) as a 
function of time, as illustrated in Fig. []Ji,c ("typical" 
meaning "away from the critical case J/A"). At short 
time, independently of J and 5, quantum and classical 

1 The difference with the classical rescaling reflects the fact that 
the norm of a quantum spin \J§? = S \J (1 + ^) is not exactly S. 
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contrast coincide that is to say C(t) « 1 — P/U with 
ti = We will therefore barely discuss this short time 
regime from now on. On longer time scale, the quantum 
contrast C(t) exhibits fast oscillations that are modulated 
by an envelope. As we are going to show, the fast oscil- 
lations essentially encode the dynamics of the single spin 
direction ni(t), or more properly the dynamics of nf(i). 
To leading order, this dynamics strongly ressembles the 
classical dynamics and is characterized by an oscillation 
time scale t Q of the order of the classical period t p (see 
Fig. [l^,c inset). By contrast, the envelope encodes the 
quantum dynamics of the single spin contrast C\(t) that 
has no classical counterpart. For large value of 5, this 
envelope exhibits a very rich multi-scales dynamics with 
many S dependent quantum time scales. In this work, 
as indicated by the arrows in Fig. [2^i,c, we concentrate 
on two of these quantum time scales: the envelope time 
t e where the envelope first vanishes and the approximate 
reccurence time t ar where it almost recovers its maximal 
value. 

Quantum spins 1/2. — The case of two spins S = \ 
can be solved exactly as the eigenvalues and eigenvectors 
are easily obtained analytically n\. The contrast is 



C(t) 



, Jt . ujt J Jt out 

cos( 7! )cos( 7! ) + - sm( 7! )sin( 7! ) 



(i) 



where w = \/J 2 + 3 A 2 . The single spin quantities C\(t) 
and Hi (t) can also be computed analytically 7 . The prod- 
uct Ci(t)fii(t) — is plotted in Fig. 3 and will be dis- 
cussed in a later part. The contrast C(t)\ given by Eq. 
exhibits fast oscillations modulated by an envelope. Be- 
cause of the latter, the contrast always vanishes at some 
time whatever the ratio J/ A. This means that there is no 
synchronization transition in this S = \ quantum case. 
However, there are still two distinct regimes separated by 
a crossover around J = A: 
(i) For J«A, the contrast is 



C(t) 



cos (At) cos( 



Jt , 

71' 



(2) 



It features fast oscillations leading to a periodic vanish- 
ing of the contrast with a period t Q = J. These fast 
oscillations correspond to a single spin direction identi- 
cal to the classical result H\(t) « [cos ( At), sin (At), 0]. 
The contrast envelope is given by the single spin contrast 
C\(t) = | cos ( 7^=) | and is characterized by an envelope 

rr- and a reccurence time t nr ~ 2t P . For a 



time t P 



2J 



generic value of J/ A the relevant Bohr frequencies are in- 
commensurate and the contrast is only quasi-periodic, 
(ii) For J ^> A, the contrast is 



C(t) 



V3A 2 3A 2 jl 
cob(— t)[l-^ani (^) 



(3) 



The fast oscillations of period t Q = ^y^, have now 
a small amplitude of order A 2 /J 2 that do not lead 
to a vanishing contrast. These oscillations are well 
in correspondance with the component |nf(t)| of the 
single spin direction Hi(i) = sign(cos(^|j-t))[l — 
|^sin 2 (^|) ) ^sin(^|)cos(^|) ) ^sin 2 (^)] (the lat- 
ter sharing the same structure as the classical result Hi 
except for the global sign function) . The contrast envelope 
corresponds to a single spin contrast C\(i) = | cos(^j-t)\ 
with an envelope time t e 



and a reccurence time 



is similar to that of off-resonance Rabi 



Vsa 

t ar ~ 2t e . The contrast still vanishes but on the much 
longer timescale t e , therefore mimicking synchronization 
during the short time dynamics. The physics behind the 
behavior t e 
flopping 6 . 

(iii) The crossover from one regime to the other occurs in 
the vicinity of J = A where the contrast has a simple ex- 
pression C(i) = |cos(^t)| 3 such that all time scales are 

of the same order t Q ~ t e 



i_ 

A ' 



Numerics. — In order to study the influence of the 
spin size, for S between 1 and 20 and various J/ A val- 
ues, we numerically obtain the eigenvalues and eigenvec- 
tors of the time-independent Hamiltonian of size (2S + 
l) 2 x (25 + l) 2 and then compute the single spin time 
evolution (5i(t)), that contains all the necessary informa- 
tion. As illustrated in Figs. []Ji,c, for S = 10, we obtain 
two regimes of behavior for the contrast C(t) depending 
on J being larger or smaller than A. For J <C A, the 
contrast has fast oscillations leading to a periodically van- 
ishing contrast with the oscillation period t Q ~ 1/A (Fig. 
[]Ji). These oscillations are modulated by an envelope with 
an envelope time that can be fitted as t e ~ y/S/J. For 
J ^> A, the contrast features fast oscillations of small am- 
plitude and with oscillation period t Q ~ 1/J (Fig. [1J). In 
this regime the contrast only vanishes with the envelope 
on a much longer time scale t e ~ y/SJ/A 2 . Extrapolating 
these two asymptotic behaviors, the ratio t e /t Q is minimal 
for J ~ A (t is maximal and t e is minimal). In fact, in 
this crossover regime J ~ A and for the studied values 
of 5, it appears difficult to correctly distinguish the two 
time scales t Q and t e . As illustrated in Figs. [IJd, after an 
initial short time decay, the contrast fluctuates a lot and 
we mainly observe that the average amplitude of these 
fluctuations decreases with increasing S. 

Large S approach. — To study whether quantum 
fluctuations destroy the classical synchronization for any 
finite S (even arbitrary large), we have developed a large 
S cumulant approach following 9 . We refer to 7 for the 
presentation of these lengthy calculations. Here, we just 
plot the contrast obtained using the cumulant approach 
and compare it with the classical dynamics and the quan- 
tum numerics in Fig. [l] (full red line). We observe that 
the cumulant approach deviates from the classical dynam- 
ics and correctly follows the quantum dynamics up to a 
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critical time t c at which it deviates from the latter by over- 
estimating quantum corrections. In the dephased regime 
J <C A we find t c ~ t e whereas in the synchronized regime 
J > A, t c ~ t e /20 . 
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Fig. 2: Contrast C as a function of rescaled time t [units of 
^] for 5 = 10. Comparison between small J s /A s (dashed 
green line), large J s /A s (full red line) approaches and quan- 
tum numerics (dotted blue line) in 3 regimes: (a) J/A = 0.2 
(dephased); (b) J/A = 1 (critical) and (c) J/A = 10 (synchro- 
nized). Initial decay (U), fast oscillations (t ), envelope (t e ) 
and approximate recurrence (t ar ) times are indicated. 

Effective Hamiltonian at small J s /A s . — In the 
limit J s /A s <C 1, the dynamics in the complete Hilbert 
space appears to be well described by the effective Hamil- 
tonian H Js ^ As = A s (5f - 5f) + J s 5f5f 7 . Its eigen- 
vectors and eigenvalues are |5i = 5, mi; 62 = 5, 777,2) and 
4i,m 2 = A s (mi - m 2 ) + J s mim 2 with = -5, . . . , 5. 
In that situation we find that the contrast is 



C(i) 



cos(Ai) cos 2S ( 



Jt 



2^5(5 + 1) 



(4) 



which recovers Eq. ([2| when 5 = 1/2. Equation Q fea- 
tures fast oscillations leading to a periodic vanishing of 
the contrast with a period t Q = -g- . These oscillations cor- 
respond to a single spin direction identical to the classical 



and spin 1/2 results: ni{t) « [cos (At), sin(At), 0]. Eq. Q 
reveals an envelope corresponding to a single spin contrast 

)S 25 (; 



Ci(t) 



Jt 



such that for large but finite 



2^S 
J 



and an en- 



*2y/S{S+l) 

5 ^> 1 it predicts an envelope time t e ~ 

velope reccurence time £ ar = 2n V s j s + }l ^ ttV St e . Com- 
parison to numerics is excellent for times t <t ar , see Fig. 
[2^i. Numerics show however that the reccurence is only ap- 
proximate, implying the existence of longer characteristic 
quantum time scales not accessible within this analytical 
approach. Note that when 5 — >• 00, then t e — >• 00 and the 
classical result C\{t) = 1 is recovered. 

Effective Hamiltonian at large J s /A s . — In the 
limit J s /A s ^> 1, using the Bloch-Horowitz projec- 
tion method 10 , we obtain an effective Hamiltonian 

Hj 3>> ^ = f [5? -25(5+1)1] + Ja2 ^ 4 V 1} [(2S)H- (S? ) 2 ] 

where St = Si + 62 • It describes the dynamics in a trun- 
cated Hilbert space restricted to the largest Bloch sphere 
of radius 2S 7 . Its eigenvectors and eigenvalues are 

,M 2 + ct with M = 



\S t = 25; M) and E M = -j^s&s-iy 
—25, . . . , 25. This approach automatically integrates over 
the fast oscillations and only captures the envelope giving 



a contrast C(t) = C\(t) = 



cos 



AS- 



V 2JS (4S-1) ) 



Using 



2JS(4S- 

a slightly different approach, which consists in keeping the 
two largest Bloch sphere (25 and 25 — 1) in a truncated 
Hilbert space 7 , we find that the contrast is actually 



C(t) 



cos 



4S-1 



AHJl 



2J(45- 1) 



) x 



[1 



A 2 (l + |) 
2J 2 



Jt 



sin 



1 



(5) 



which recovers Eq. ([3| when 5 



tures fast oscillations with period t Q 

A 2 (l+-1) 



1/2. Equation ^ fea- 

TTy/l + l/S 
~ J 



2J 2 



and small 
that coincide with the clas- 



amplitude of order 
sical result up to 1/5 corrections. For large but finite 
5 ^> 1, the single spin contrast C\{t) predicts an enve- 
lope time t e « 4 



2ttJS(4S-1) 



25 7 and an envelope reccurence time 
2St e . As in the previous case, the 



A 2 
- 7T 



agreement with numerics is very good for times t < i ari 
see Fig. [2J3. More importantly, the main result of this 
section is that the contrast first vanishes at the envelope 
time t e ~ >/5J/A 2 , that diverges when 5 — >> 00, there- 
fore pointing toward synchronization in the classical limit. 
Nevertheless, for large but finite 5 ^> 1, the quantum as- 
pect is striking when considering the norm of the total 

spin: on the one hand \J (5|) = ^/2S(2S + 1) is a con- 
stant, as the system remains on the largest Bloch sphere, 
but on the other \(S t )\ goes to zero at the envelope time. 
These quantities correspond to two definitions of the norm 
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of a vector which coincide classically but not quantum me- 
chanically. 

As illustrated in Fig. [2}}, we note that neither Eq. Q 
nor Eq. ([5| correctly captures the quantum contrast 
around the critical point J/ A = 1 of the classical dy- 
namics. 




Fig. 3: Dynamics of a single spin (Si) / S = Ci(t)n\(i) on the 
Bloch sphere for times < t < t ar - Comparison between quan- 
tum (dotted blue line) and classical (full red line) dynamics in 
3 regimes: First column (a,d): J/A = 0.2 (dephased); second 
column (b,e): J/ A — 1 (critical); third column (c,f): J/A = 10 
(synchronized). First line (a,b,c): S = 10; second line (d,e,f): 
S = 1/2. The equator is indicated with a black line. 

Single spin on the Bloch sphere. — The single 
spin representation (Si)/S = Ci(t)ni(t) is particularly 
powerful in providing a physical picture for the system's 
dynamics and in comparing the classical and the quan- 
tum cases. Figure [3] show the single spin behavior on the 
Bloch sphere in three regimes - dephased J <C A, critical 
J = A and synchronized J ^> A - for both the classical 
(red lines) and quantum dynamics (blue dots) and for two 
spin sizes S = 10 and S = 1/2. Away from J/A = 1, the 
dynamics of the single spin direction ni(t) is quite simi- 
lar in the classical and quantum cases: in the dephased 
regime, it stays almost in the equatorial plane and turns 
around z axis with a period t Q ~ tt/A (see Fig. [3^i,d); in 
the synchronized regime, apart from small (S dependent) 
amplitude oscillations of period t Q ~ tt/J, the direction 
n\(t) is almost locked and aligned with +x except that 
in the quantum case it jumps back and forth from +x to 
— x at every odd multiple of the approximate reccurence 
time t ar (see Fig. [3)3, e and the sign function in n\ below 
Eq. ([3])). Most of the quantum behavior is thus encoded 
in the dynamics of the single spin contrast C\(t) which 
vanishes at the envelope time t e and almost recovers at 
the reccurence time t ar , where the S dependent timescales 
t e and t ar depend on the considered regime. In the de- 
phased regime the combined dynamics of C±(t) and n\(t) 



results in an inwards spiral motion for the single spin in 
the equatorial plane (see Fig. [3^i,d) until t e , followed by 
an outwards spiraling until t ar . For S = 10, the long time 
interval t e < t < t ar during which C\(t) ~ corresponds 
to the higher density of points visible in the center of the 
Bloch sphere in Fig. [3^i. In the synchronized regime, the 
combined dynamics now results in a cigar-like shape of 
the complete trajectory visible on Fig. [3j3,f. In the critical 
case, the classical direction n\ goes monotonically from x 
to z. The quantum dynamics is here quite different and 
depends strongly on S. More precisely, when increasing 
5, it evolves towards the classical one for times shorter 
than the envelope time (see Fig. |3]3,e), however at longer 
times, the quantum single spin always ends up visiting the 
inner part of the Bloch sphere, while the classical spin is 
constrained to its surface. 

Qualitative picture. — The main qualitative picture 
emerging from our study is that the dynamics of a quan- 
tum spin (and of the quantum constrast) results from the 
combined dynamics of its effective direction n\ (t) and its 
effective norm or envelope C\(t). On the one hand, the 
short and fast time dynamics are governed by the direc- 
tion n\(t) and, apart from 1/5 corrections, it retains most 
of the classical aspect of the spin dynamics; in particular 
ti = ^ is the same in the quantum and classical cases 
and does not depend on J or S. In addition, away from 
J = A, the fast oscillations time t Q ~ min( J, j) cor- 
responds well to the classical period t p . On the other 
hand, the slow and long times quantum spins dynamics 
are governed by the envelope C\ (t) . Since C\ (t) is a direct 
measure of the effective spreading of the spin wavepacket 
on the Bloch sphere, this part of the dynamics is a pure 
quantum phenomenum and it strongly depends on S. For 
large value of S, this envelope C\(t) has a very rich mul- 
tiscale dynamics featuring collapses (C\ ~ 0), revivals 
and reccurences (C\ ~ 1). Quantitatively, we found that, 
away from J = A, the envelope time of the first collapse 
t e ~ y/S max( j, -j^) and the time of the first approximate 
recurrence t ar ~ S max( j, ^2) ~ \/St e both diverge with 
S in the classical limit. More generally, we expect an 
increasing hierarchy of collapse and reccurence quantum 
timescales with increasing values of S. 

The way the classical limit emerges from the quantum 
dynamics as S —> 00 is thus quite remarkable. The clas- 
sical dynamics is characterized by only two timescales (ti 
and t p ~ t ), whereas the quantum dynamics has more 
and more timescales as S increases (t e , i ari etc.), so that 
a priori it seems unlikely that the two pictures would rec- 
oncile in the large S limit. However all of the quantum 
timescales diverge when S —> 00 so that in the end, the 
classical dynamics emerges as the short time behavior of 
the quantum dynamics. The transition occurs between 
the oscillation to and the envelope t e times. The situation 
close to the classical critical point J/A = 1 is quite fa- 
vorable to observe quantum effects, as the classical period 
diverges (see Fig. |2b) and the envelope time is minimal, 
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although it is difficult to study analytically. 

A similar behavior - namely multi-scale quantum dy- 
namics and phase transition in the classical dynamics 
only - was found in a simpler non-linear single large spin 
model. This is the Bose-Hubbard dimer, describing e.g. a 
Bose-Einstein condensate in a double-well potential 12 



and which, in the spin language, maps onto the Lipkin- 
Meshkov-Glick model, see e.g. [13]. In the classical or 
mean-field limit the condensate is described by a Gross- 
Pitaevskii equation, which features a self-trapping transi- 
tion LL4] . At small onsite repulsion, the bosons initially all 
in the same well can tunnel to the other well. But when 
the onsite repulsion exceeds a critical value, the bosons are 
essentially trapped in their initial well. Including quantum 
fluctuations has a dramatic effect: the self-trapping tran- 
sition no more exists - for all interaction strength, the 
condensate can visit both wells - and a rich multi-scale 
quantum dynamics emerges, see e.g. [15] . The complete 
dynamics involves at least three very different timescales, 
the two largest being quantum and increasing very fast 
with the spin size [15] . 

Experimental relevance. — In 0, a spin kinetic 
equation in energy space S(E,t) was developed and solved 
numerically in order to compare to the experimental re- 
sults. The agreement between theory and experiment 
was quite good but the spins were intrinsically classical. 
Here, the two spins model captures the physics of the 
synchronization transition in the large S limit and also 
includes the effect of quantum fluctuations. In the de- 
phased regime, it features spurious large amplitude fast 
oscillations, which are not observed in experiment and are 
due to the huge simplification in the size of the Hilbert 
space from the microscopic dimension 2^ to the effective 
dimension (25 + l) 2 = (TV + l) 2 . The most interesting as- 
pect of the two spins model is the appearance of quantum 
dynamics on the timescale of the envelope t e . From the 
experiment R , we get A ~ 10 5 rad/s, J ~ 10 5 rad/s and 
S = N/4 ~ 10 4 so that the envelope time t e ~ t e S ~ 10 s 
is comparable to the duration of the experiment (~ 5 s). It 
would be interesting to probe longer timescales experimen- 
tally to see the appearance of quantum dynamics modu- 
lating the classical behavior by sequences of collapses and 
revivals. The two spins model is not only relevant to the 
atomic clocks |4,5 1, but should also concern other experi- 
mental systems such as two quantum dots trapping spins 
[16] , two nanodisks with large magnetic moments [17] or 
generally systems with two coupled large spins such as 
molecular magnets, in particular so-called molecular mul- 
tidot devices [18] . 

Conclusion and perspectives. — In conclusion, we 
introduced and studied a quantum model of two exchange- 
coupled spins S in the presence of an inhomogeneous mag- 
netic field. Our main result is that quantum fluctuations 
prevent the occurrence of the synchronization transition 
found in the classical limit. The quantum dynamics is 



very rich, featuring collapses, revivals and recurrences de- 
scribed by multiple timescales. The most important quan- 
tum feature is the envelope of the contrast, which corre- 
sponds to the spreading of a single spin wavepacket. The 
classical dynamics is reached in a non-trivial way emerg- 
ing as the short time behavior of the quantum dynam- 
ics. The quantum dynamics at the classical critical point 
J = A is especially intriguing and deserves further inves- 
tigations. In the future, we plan to study the possibility of 
spin squeezing 19 within our simple model. It would also 
be interesting to include the effect of coupling to a bath 
that could favor a synchronization transition by rendering 
the system more classical [20], 
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